Non-adiabatic level crossing in (non-) resonant neutrino oscillations 
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We study neutrino oscillations and the level-crossing probability Plsz = exp(— ■y n J r n ir/2) in 
power- law like potential profiles A(r) oc r n . After showing that the resonance point coincides 
only for a linear profile with the point of maximal violation of adiabaticity, we point out that 
the "adiabaticity" parameter j n can be calculated at an arbitrary point if the correction function 
T n is rescaled appropriately. We present a new representation for the level-crossing probability, 
Plsz = exp( — K n Q n ), which allows a simple numerical evaluation of Plsz in both the resonant and 
non-resonant cases and where Q„ contains the full dependence of Plsz on the mixing angle #. As 
an application we consider the case n — —3 important for oscillations of supernova neutrinos. 
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Introduction — The analytical study of non-adiabatic 
neutrino oscillations has a long history. Soon after the 
discovery of "resonant" neutrino flavor conversions by 
Mikheev and Smirnov [Q, the leading non-adiabatic ef- 
fects were calculated for a linear potential profile in 
Rcf . 0] . They were obtained in the form of the Landau- 
Stiickelberg-Zener (LSZ) crossing probability M 



-Plsz = exp 



The adiabaticity parameter 7 for a linear profile is 
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where E is the energy, A = m 2 , — m\ > 0, S = sin(2$), 
C = cos(2$), 1} G [0 : 7r/2] is the (vacuum) mixing an- 
gle and A = 2EV = 2y/2G F N e E is the induced mass 
squared term for the electron neutrino. The parameter 
7 has to be evaluated at the so-called resonance point 
ro, i.e. the point where the mixing angle in matter is 
i9 m = 7r/4. For a linear profile, adiabaticity is maximally 
violated for i) m — 7r/4 and, therefore, the probability 
that a neutrino jumps from one branch of the dispersion 
relation to the other is indeed maximal at r Q . 

Later, Kuo and Pantaleone derived in Ref. j| the LSZ 
crossing probability for an arbitrary power-law like pro- 
file, A oc r n . They found that also in this case the de- 
pendence on the neutrino masses and energies can be 
factored out, while the effect of a non- linear profile can 
be encoded into a correction function J- n , 



Plsz = exp (-^ T n {d) 



(3) 



which only depends on •& and n. The "adiabaticity" pa- 
rameter 7„ has to be evaluated still at ro although, as 
we will show, it does not coincide with the point of max- 
imal violation of adiabaticity (pmva) for n 7^ 1. An un- 
satisfactory feature of Eq. (||) is its restricted range of 



applicability: the resonance condition i9 m = 7r/4 has a 
solution only if d < ir/4 for neutrinos or if -d > ir/4 for 
anti-neutrinos, respectively. Therefore, it is not possi- 
ble to calculate analytically, e.g., the survival probability 
of supernova anti-neutrinos in the quasi-adiabatic regime 
assuming a normal mass hierarchy. 

The purpose of this Letter is twofold. First, we dis- 
cuss the physical significance of the resonance point ro- 
We show that the product 7„J 7 „ can be evaluated at an 
arbitrary point and conclude that the "resonance" point 
ro has, for a general profile n 7^ 1, no particular physical 
meaning: neither does it describe the point of maximal 
violation of adiabaticity nor is it necessary to calculate 
the "adiabaticity" parameter j n at r Q . Second, we pro- 
pose a new representation for Plsz that is valid for all ■& 
and allows an easy numerical evaluation. As an applica- 
tion, we consider the case n — — 3 which is important for 
oscillations of supernova neutrinos. 

Resonance point versus point of maximal violation of 
adiabaticity — We use as starting point for our discussion 
the evolution equation for the medium states ip first given 
in Ref. @, 



d_ 

d7 

Here. 



t 1 



iA m /(4£) 



AA B %E))(f 2 ) Nl 



(A - AC*) 2 + (AS) 2 



(5) 



denotes the difference between the effective mass of the 
two neutrino states in matter, $ m is the mixing angle in 
matter with 



tan 2§„ 



AS 



AC- A 



(6) 



and d' m = d-d m / 'dr. 

The evolution of the neutrino state is adiabatic at a 
given r, if the diagonal terms are large compared to the 



1 



off-diagonal ones, |A m | ^> |4i?z9J n |. Thus the point where 
adiabaticity is maximally violated is given by the mini- 
mum of A m /i}' m . Differentiating H 



1 



2A 2 S 2 



sin 3 2tf m dA/dr 
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for the case of a power-law profile, A(r) oc r n , we find 
the minimum at 



cot(2i? m - 

- [cot(2i?. 
n 



2i9) + 2cot(2tf m ) 
n - 21?) - cat(20„)] 



0. 



(8) 



For n = l, the pmva is at i? m = 7r/4 for all Thus, 
in the region where the resonance point is well-defined, 
the pmva and ro coincide. In the general case, n / 1, 
the pmva however agrees with the resonance point only 
for $ = 0. Finally, we recover the result of Ref. || for an 
exponential profile in the limit n — > ±oo. 
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FIG. 1. Change of the survival probability dp(r)/dr of 
a neutrino produced at r = as V2 together with the point 
of maximal violation of adiabaticity (dots) and the resonance 
point (stars) for a power-law profile A oc r -3 . The height of 
the different curves is rescaled. 

In Fig. ^, we show the change of the survival proba- 
bility dp(r)/dr = d|"02(?*)| 2 /dr of a neutrino produced 
at r = as i>i as function of r together with the point 
of maximal violation of adiabaticity predicted by Eq. (0) 
and the resonance point for a power-law profile A cx r . 
It can be clearly seen that Eq. (JsJ) accurately describes 
the most probable position of the level crossing, while 
the resonance point predicts a transition in less and less 
dense regions until for $ — it/ 4 the concept of a resonant 
transition breaks down completely. If the true profile is 
only approximately a power-law, its exponent should be 
determined therefore by the region around the pmva, not 
by the region around the resonance point. 



The correction functions T n — We now recall briefly the 
calculation of the leading term to the crossing probabil- 
ity within the WKB formalism || . In the ultrarelativistic 
limit and omitting an overall phase, the WKB formula 
results in 



IiiPlsz = --=9 



xi(Ai) 



da; A n 



(9) 



where A2 = Ae 2 " 5 is the branch point of A m in the up- 
per complex x plane. We identify the physical coordi- 
nate r G [0 : 00] with the positive part of the real x axis, 
i.e. we consider a neutrino state produced at small but 
positive x propagating to 1 = 00. Then, a convenient 
choice for A\ is to use the real part of A2 for C > 0, 
i.e. the "resonance" point Ai = AC. However, we stress 
that this choice has technical reasons and makes sense 
only for C > 0: consider for instance the simplest case 
n = 1. Then both the integration path chosen and the 
branch cut are for C < in the half-plane $t(x) < 0. The 
physical interpretation is therefore that an anti-neutrino 
state created at small but negative x propagates to —00. 
This case is however equivalent to a neutrino state prop- 
agating with C > in the right half-plane and therefore 
contains no new information. Thus, we expect the cor- 
rection functions J- n obtained with the integration path 
from AC to Ae to be functions with period 7r/4 and 
to be valid only in the resonant region. 
After the substitution A — A(B + C), 
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is 



dB — (B 2 + S 2 ) 1 / 2 , (10) 
dB 



one has to expand the Jacobian dx/dB in a power series 
in order to solve the B integrals. Kuo and Pantaleone 
chose as expansion point the "resonance" point Bq = 0, 
because it leads to the simplest result. It is this choice, 
arbitrary from a physical point of view, that leads to the 
evaluation of the parameter 7„ at Bq = 0. In general, a 
change of Bq in the definition of 7„ will be compensated 
by such a change in the correction function T n that the 
physical observable Plsz is independent of Bq . 

The final result for the correction functions given by 
Kuo and Pantaleone was 
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This series for JF„(i?, 0) converges only for ■& < tt/8 and is 
therefore not suited, even in the resonant region, for nu- 
merical evaluation in the phenomenologically most inter- 
esting case of maximal or nearly-maximal mixing. Rep- 
resenting the series as a hypergeometric function, 



^(0,0) = 2 Fi 
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tan 2 (2tf) , (12) 
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one can use, however, the Euler integral representation g] 
of 2-P1 as the analytical continuation for all i) £ [0 : 7r/4]. 

We now discuss the non-resonant case, ^ £ [f : §]• 
The expression 7„(0).P n (0, $) becomes ill-defined in this 
region for several reasons: the assumption C > used in 
the derivation for J r n ('d ) 0) is not fulfilled, and the "reso- 
nance" condition B = or A = AC has no solution for 
i? € [7r/4 : tt/2}. Moreover, the integration path used in 
Eq. (10) can be for C < in the unphysical region x < 0. 
However, we know from the results above that the cross- 
ing probability is non-zero also in the non-resonant case. 
Also, the WKB method should work independently of the 
sign of C as long as the evolution of the neutrino state is 
not strongly non-adiabatic. 

We now show that it is not necessary to evaluate i n Tn 
at the resonance point Bq = 0. Requiring the invariance 
of -Plsz against variations of Bq, 



Tn 



din 

dB 



In 



dB Q 



0. 



we obtain with 



dB n(P + C) 
a differential equation for AT /&B$. Its solution 



c 



B + C 



l/r, 



(13) 



(14) 



(15) 



allows us to rescale the correction functions obtained for 
Bq = to arbitrary Bo. Therefore, we can calculate Plsz 
now also for the non-resonant region in the two cases in 
which the function T n is known for all 

For a l/r profile, A(r) = 2EVa(Ra/r), the correction 
function is jl| 



(l-tan 2 (tf)) ; 
1 + tan 2 (tf) 



(16) 



If we want to calculate the level-crossing probability in 
both the resonant and non-resonant cases, we have to 
choose Bq > 1, e.g. 



^(1,^ = 1^^(0,^); 
the crossing probability follows as 

Plsz = exp {-2n R V sin 2 (z9)} 



(17) 



(18) 



As an example, we compare in Fig. g the results of a 
numerical solution of the Schrodinger equation (^) with 
the analytical calculation of the Plsz using the rescaled 
F_! function ||. The agreement between the different 
methods is excellent. 

In the limit n — > ±oo, which corresponds to an expo- 
nential potential profile, the scale factor [C / (B + C)] 1 /™ 



goes to 1 and the correction function becomes indepen- 
dent of Po: as it should be. The resulting crossing prob- 
ability is 



P 



LSZ 



exp 



vrAPo 
E~ 



sin 2 (tf) 



(19) 



In Rcfs. 0, this expression was derived by solving 
Schrodinger's equation directly. In these works, it was 
assumed that the obtained expression is valid only in the 
resonant region and only recently was it pointed out that 
it is valid for all $ [0. Note also that In Plsz has the 
same, very simple dependence on for both n = — 1 and 
n = ±oo. 
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FIG. 2. Crossing probability P c for A oc l/r with 
RoVo = 0.2: numerically (solid line), with T-\ evaluated at 
Bo = 1 (squares) and Q-\ (circles). 

The correction functions Q n — We start directly from 
Eq. (^|), but use now as integration path in the com- 
plex x plane the part of a circle centered at zero starting 
at Ai — A and going to the end of the branch cut at 



A 2 = Ae AW . Substituting x = R (A/A ) 1/n e^ 



the 



case of a potential A = Ao(r / Rq) 11 , we can factor out the 
i? dependence of Plsz into functions Q ni 
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where 



A 



l/n 



is independent of i? and 

r 2i?/r, 



/ dip e iv (e mip - C) + S' 2 
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(20) 



(21) 



(22) 



The functions Q n are well suited for numerical evaluation 
and always correspond to a neutrino state propagating in 



3 



the physical part of the x plane, x > 0. Therefore, they 
have, in contrast to the T n functions, the period tt/2 and 
are valid for all i?. In Fig. ||, the results of this new rep- 
resentation are compared with those obtained above for 
the case n = — 1. 

Oscillations of supernova neutrinos — The potential 
profile A(r) in supernova (SN) envelopes can be ap- 
proximated by a power law with n w —3, and V(r) — 
1.5 x 1(T 9 eV (10 9 cm/rT 3 |l2). Since only v e were 
detected from SN 1987A and also in the case of a fu- 
ture galactic SN the v e flux will dominate the observed 
neutrino signal, an analytical expression for P c valid in 
the non-resonant part of the mixing space is especially 
useful [0. The probability for a v e to arrive at the sur- 
face of the Earth can be written as an incoherent sum of 
probabilities 

P gg = (l-P c )cos 2 i? + P c sin 2 tf. (23) 



an "adiabaticity" parameter 7„ evaluated at the reso- 
nance point and a correction function J- n is misleading 
for n ^ 1: neither does the level-crossing probability 
have a maximum at ro nor does this splitting allow the 
calculation of Plsz in the non-resonant region. We have 
proposed a new representation for Plsz that avoids these 
problems and is hopefully useful for the investigations of 
oscillations of supernova neutrinos. 
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FIG. 3. Contours of constant survival probability Pes, 
numerically (solid lines), with Q-z (dashed lines) and J--z 
(squares, only for i? > 7r/4), for A tx r~ 3 as given in the text. 

In Fig. |^, we compare the results of a numerical solu- 
tion of the Schrodinger equation (|j) with the analytical 
calculation of P gs using the £/_3 and the T--$ function. 
The latter is shown only for its range of applicability, 
d > 7r/4. The agreement between the two methods us- 
ing the WKB approach is again (for > 7r/4) excellent. 
Generally, these two methods agree also very well with 
the results of the numerical solution of the Schrodinger 
equation; there are only small deviations in the regions 
where the contours change their slope. 

Summary — We have discussed non-adiabatic neutrino 
oscillations in general power-law potentials Aar". We 
have found that the conventional splitting of In Plsz in 
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